Thermodynamic limits of dynamic cooling 
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We study dynamic cooling, where an externally driven two-level system is cooled via reservoir, a 
quantum system with initial canonical equilibrium state. We obtain explicitly the minimal possible 
temperature Tmin > reachable for the two-level system. The minimization goes over all unitary 
dynamic processes operating on the system and reservoir, and over the reservoir energy spectrum. 
The minimal work needed to reach Tmin grows as 1/Tmin- This work cost can be significantly reduced, 
though, if one is satisfied by temperatures slightly above Tmin- Our results on Tmin > prove 
unattainability of the absolute zero temperature without ambiguities that surround its derivation 
from the entropic version of the third law. The unattainability can be recovered, albeit via a different 
mechanism, for cooling by a reservoir with an initially microcanonic state. We also study cooling via 
a reservoir consisting of A*' 3> 1 identical spins. Here we show that Tmin ex: i and find the maximal 
cooling compatible with the minimal work determined by the free energy. 

PACS numbers: 05.30.-d, 05.70.-a, 07.20.Mc 



I. INTRODUCTION 

Some physical systems have to be cooled before they 
can demonstrate interesting features, e.g. quantum prop- 
erties of matter are typically displayed only after suit- 
able cooling. There are various cooling methods; we dis- 
tinguish here the "brute force method", the Nernst-set 
up, and the so-called dynamic cooling. The brute force 
method of cooling amounts to bringing the system in 
contact with a low-temperature thermal bath, so that it 
relaxes to this lower temperature. 

When trying to do without a pre-existing cold bath 
different strategies are needed: Low temperatures can 
alternatively be produced dynamically from an initially 
equilibrium system via cyclic action of an external field 
|ll-[l8J|. In this case the resulting cooled state should be 
exploited before the system has time to relax back to 
equilibrium. 

Within macroscopic quasi-equilibrium thermodynam- 
ics some of these methods have been summarized in the 
Nernst set-up of cooling, which supports a formulation of 
the third law X'S]- In fact, this law controls the cooling 
of a macroscopic target initially in contact with an ap- 
propriate macroscopic reservoir. The target is then typ- 
ically subject to a two-step cooling process, in which an 
external field cycle is executed under (quasi-equilibrium) 
isothermal and adiabatic conditions, respectively. A well- 
known realization of this Nernst set-up is the magne- 
tocaloric effect which was first observed in 1880, but still 
attracts attention [8|. 

There are however other methods of cooling that be- 
came important with the rise of low-temperature physics 
|9l- [l3 . [l4| - [l8| . Here the target of cooling (or the reservoir, 
or both) is generally not a macrosocpic system, while the 
process that produces low temperatures ceases to be a 
quasi-equilibrium one. Hence the understanding of these 
methods should rely on the actual dynamics rather than 
on quasi-equilibrium thermodynamics. We therefore re- 



fer to them as dynamic cooling. A notable example of 
this is the dynamical nuclear polarization in NMR [2- 
[la . [l4| . Within this method the nuclear spins are cooled 
via transferring polarization from electron spins by means 
of external microwave fields J9l-[l2l. [l4| . It was originally 
employed in the solid-state NMR, but since recently it 
is also applied for the liquid-state NMR in view of its 
medical and bio-physical applications [Ij]. Other exam- 
ples of dynamic cooling are algorithmic cooling [1^ [l^ , 
bath-assisted cooling j[17i] and laser cooling of motional 
state in atoms |18| . 

Dynamic nuclear polarization illustrates the basic in- 
gredients of dynamic cooling processes that are seen al- 
ready in the Nernst set-up: the target system to be cooled 
(nuclear spins), the reservoir, which plays the role of a 
polarization source or the entropy sink (electrons), and 
external fields that couple these two together (microwave 
radiation at suitable frequency) . Recall that no cooling is 
possible without reservoir [l9|; i.e., one can never cool the 
whole system {— target -I- reservoir) by coherent fields. 

As already indicated, the operation of cooling and 
the concept of the third law are intimately interrelated. 
While we do not attempt to dwell on the general ideas be- 
hind operationalism as originally promoted by Bridgman 
|20| . it is nevertheless worthwhile to stress the potential 
benefits of such an approach even to modern thermody- 
namics, and, in particular, for reaching the objectives of 
this paper. Our aim here is threefold: To formalize cool- 
ing schemes in terms of imposed limited resources, to 
parametrize specific reservoir models, and to investigate 
the resulting minimal temperature as a (scaling) function 
of those operational parameters. 

In particular, we intend to study the minimal temper- 
ature Tmin reachable within a sufficiently general set-up 
of dynamic cooling, where both the target and reservoir 
arc finite quantum systems, and where the final state of 
the reservoir can generally be far from equilibrium. We 
determine how Tmin depends on resources of the set-up 



and note that the minimal work necessary to achieve Tmin 
grows as yi— whenever Tmin — ?► 0. Hence the work is a 
diverging resource of dynamic coohng, in contrast to the 
Nernst set-up and the third law, where the work done for 
cooling is always well-bounded. We also determine the 
minimal temperature (and the minimal work needed to 
attain for it) for several concrete reservoir models. 

Our results on T^-m > extend (to dynamic cooling) 
the unattainability formulation of the third law. Hence 
section in] reviews the (frequently disguised) assumptions 
of the Nernst set-up and of the third law in its two for- 
mulations. (The reader who is well aware of the third law 
may just consult section [HP I for relevant assumptions.) 

The dynamic cooling set-up is defined in section IHII 
Section FHI CI shows that dynamic cooling with a reservoir 
in an initial canonical state does not allow to reach the 
absolute zero of temperature [unattainability] . In section 
IHIEI we summarize relevant conditions of the dynamic 
cooling set-up. 

The following three sections are concerned with reser- 
voir models. Section ITVl studies the maximal cooling pos- 
sible within the dynamic set-up. It confirms that the 
work is a relevant resource for cooling and determines the 
minimal work necessary to achieve the maximal cooling. 
Section|V]studies a reservoir with the homogeneous, non- 
degenerate spectrum, while in section IVll the reservoir is 
modeled as a thermal bath consisting of A^ ^ 1 spins. 
This model is complemented by a discussion on the dy- 
namic cooling process for initially microcanonic state of 
the reservoir. 

In the last section we summarize our results and dis- 
cuss their relations with other microscopic approaches for 
studying thermodynamic limits of cooling: the approach 
of Ref. [2l| , where the cooling process is restricted to op- 
erations within degenerate subspaces of the joint Hamil- 
tonian of the system and non-equilibrium reservoir, and 
the approach based on quantum refrigerators [23l - l27| . We 
close by pointing out some open issues. 



II. OPERATIONAL ANALYSIS OF THE 
NERNST SET-UP AND THE THIRD LAW 

A. Description of the Nernst set-up 



We are given a macroscopic system in contact with a 
much larger thermal bath [l| . Both have initial tempera- 
ture Tin. One now switches on an external field g acting 
on the system and changes it slowly from its initial value 
f^in to its final value gtin- The contact with the bath is 
held fixed. Hence this is an isothermal process. Let gin 
and gun be chosen such that the entropy S[T,g] of the 
system decreases (see Fig. [T|) 



S[Tin,gin\ > S'[Tin,5fin] 



(1) 



The entropy difference is transferred to the bath. 

Next, the system is thermally isolated (i.e., isolated 
from the bath) and the field g is slowly returned back to 



its original value gin thereby completing the cooling cy- 
cle. (Note: the cyclic condition only refers to the exter- 
nal field, i.e. the time-dependent Hamiltonian; it would 
be meaningless to require the target to return to its ini- 
tial state.) Hence this part of the process is thermally 
isolated and reversible (adiabatic). Since the system is 
macroscopic, the adiabatic part is characterized by con- 
stant thermodynamic entropy 1]: 



S[Tf^n,gin] = S [Tin, gun], 



(2) 



where Tgn is the final temperature of the system. Due 
to H]) and 1^ > 0, we get Tfin < Tin: some cooling has 
been achieved ^^. Recall that |^ > is one of basic 
conditions of local thermodynamic stability (positivity of 
the specific heat) [1|. The system will stay cold as long 
as it is well isolated from the bath. 

This cooling cycle has a work-cost. In terms of the free 
energy 



F[T,g]^U[T,g]~TS[T,g] 



(3) 



of the system, where U[T,g] is its energy, the work done 
in the isothermal step is F[Tin, gan] — F[Tin, gin]- The 
work in the adiabatic step is U[Tf\n, gin] — C^[Tin, <?fin]- 
Summing them up and using ([2]) one has for the work 
W invested per cooling cycle ^: 



W = F[Tfin,gin] - F[Tin,gin] > 0, 



(4) 



which is always positive due to ^ — —S < and 
Tin > Tfln- Below, in sections [Hi] and EH we relate (g]) 
to the minimal work-cost of the dynamic cooling set-up. 
Note that W according to @ is always finite, even for 
Tfin = 0. 



B. Entropic formulations of the third law 

The entropic formulation states that for each fixed and 
finite value of the field g, the entropy S[T, g] > smoothly 
goes to zero with temperature T: S[T,g] — >■ for T — >■ 
[l|-|7[; see Fig. [TJ The entropic formulation is frequently 
deduced from the non-degeneracy of the system's ground 
state. However, this does not suffice for a derivation, 
because for a macroscopic system the appropriate order 
for taking limits is that first the thermodynamic limit 
is taken, the entropy density is calculated, and then af- 
ter that the low temperature limit goes [J, |5[ . Then the 
zero-temperature entropy starts to depend on global fea- 
tures of the energy spectrum. The full derivation within 
statistical mechanics was attempted several times \^j^, 
but is regarded to be an open problem, because there 



^ We use the sign convention W > 0, if the work-source looses 

energy. 
^ We assume the convention, where the entropy is defined to be 

non-negative. 




FIG. 1: A schematic representation of the relation between 
the entropic formulation of the third law and the attainability 
of r = 0. S and T denote entropy and temperature, respec- 
tively; g is the external field. Bold (dashed) lines represent a 
situation, where the entropic formulation holds (is violated). 
Arrows represent the isothermal and adiabatic parts of the 
cooling cycle. We see that the final temperature T^n is greater 
than zero. However, it is formally possible to attain T = 
(dashed arrows) when the entropic formulation is violated in 
the sense that S[T — 0] still depends on the field g. 



are different classes of theoretical Hamiltonian systems 
violating the entropic formulation: the zero-temperature 
entropy for them is not zero and depends on external 

fields 0, i H m, 113, MM ■ 

It is believed that Hamiltonians which would violate 
the entropic formulation are unstable with respect to 
small perturbations [3, l6|. Such perturbations would 
have to be generated by some appropriate noise model. 
Unfortunately, a general recipe for selecting such a model 
does not seem to be known (or even exist) ^. There are 
results showing how one should select the noise model for 
particular cases only [2, [a, IM, [4^ i4^ i49|J . 



C. Unattainability formulation of the third law 

Already the entropic formulation of the third law can 
be brought in contact with operational requirements. 
This is even more so for the unattainability formulation. 
This formulation states that it is impossible to reach ab- 
solute zero of temperature, T = 0, by finite means, i.e. 
via a tractable physical process. Note that the unattain- 
ability formulation is clearly different from the entropic 
formulation that refers to a limiting feature of a defi- 



nite function of temperature (entropy). The unattain- 
ability formulation relies on the notion of tractability, 
whose general formalization appears to be hard to come 
by. Hence the two formulations cannot be completely 
equivalent 13|. 

Within the Nernst set-up one naturally assumes that 
for a tractable physical process the strength g of the 
external field should not assume infinite values [l|, |2[. 
This suffices for a heuristic derivation of the unattain- 
ability formulation from the entropic formulation; see 
Fig. [1] and [i|. Indeed, Fig. [1] shows that also shows 
that the unattainability is violated together with counter- 
examples of the entropic formulation ^. But there are also 
examples that would violate the unattainability formula- 
tion only [3, [22| ^ . Again, one way out of this ambiguous 
situation is to look for via "sufficiently" stable Hamilto- 
nians 0]. 



D. Summary of conditions 

We finally try to summarize the Nernst set-up in terms 
of two lists of conditions: The first list Nl— N6 specifies 
details of the task: it reflects our choice for implementing 
a satisfactory cooling process. In part these conditions 
may be tied to the actual technology. The second list 
Tl— T3 introduces fundamental tractability conditions. 
These stress the unavoidable finitness of the resources 
at our disposal and should remain intact within any fu- 
ture technology. The respective process will be called 
tractable insofar as these requirements can be fulfilled. 

Nl. Initial state: equilibrium for the system and bath. 

N2. System (target of cooling): macroscopic. 

N3. Cooling process: Quasi-equilibrium. 

N4. Bath: macroscopic, much larger than the system. 

N5. External field: changes cyclically. 

N6. The Hamiltonian of the target: stable with re- 
spect to perturbations. 

Nl recalls that we intend to start from a given fi- 
nite temperature, while N2— N4 are assumed because the 
process is to be embedded into macroscopic equilibrium 
thermodynamics. N2 and N3 reflect the technology of 
Nernst's times. N5 is assumed, because the cooled sys- 
tem has to be autonomous, i.e. once it was cooled, it 



^ The reader may consult controversies with the validity of the 
entropic formulation in the field of the Casimir effect [44|446ll or 
for black holes |47| , criteria of this formulation for spin systems 
[48l |49| . and classes of one- dimensional counter-examples to the 
formulation [SOl l5l| . The authors of |5ll | ensured robustness of 
their counter-examples to a class of noise models. 



This would not imply that one can ever verify reaching T = 0, 
since any temperature measurement has a finite precision. (The 
situation with reaching a positive temperature, say 10 K, is dif- 
ferent, because one can arrange for passing through this tem- 
perature at some [uncertain] time.) In particular, temperature 
fluctuations can prevent a precise determination of low temper- 
atures l42l |43I . According to the standard thermodynamic ap- 
proach [42l | these fluctuations will grow for T — > due to the 
vanishing heat capacity. The issues of temperature fluctuations 
was recently reconsidered and clarified in [43l | from the first prin- 
ciples of quantum mechanics. 

The mechanism of these counter-examples is that the entropy 
nullifies at some Tc(g) > and stays zero for T < Tc{g). 



can be used in other places, without being kept under a 
constant external field. N6 is included here because a 
general noise model for perturbations is not known: the 
concrete stability requirement thus becomes a matter of 
choice. 

We now turn to the tractability conditions; this list is 
not necessarily complete (cf. section IVICj) . 

Tl. External field has finite strength. 

T2. The work-cost necessary for cooling is finite. 

T3. Duration of the cooling cycle is finite. 

Tl turns out to be essential for deducing the unattain- 
ability formulation of the third law. It is also relevant 
for the entropic formulation. (It is easy to find examples, 
where the entropy does not nullify with the temperature 
if simultaneously the field goes to infinity 2].) Due to ^, 
T2 is guaranteed by N1-N4. Typically, T3 becomes 
essential for certain scenarios: Indeed, in frustrated sys- 
tems reaching the true equilibrium may demand unreal- 
istically long observation times; on practically relevant 
observation times the system may find itself trapped in 
quasi-stationary states, where the non-equilibrium en- 
tropy does not go to zero with the ambient tempera- 
ture, though the true equilibrium state does satisfy the 
entropic formulation. This effect was observed experi- 
mentally |52| |. It is also well known for glassy systems 
[53] (residual entropy); see [S^ for a review. 



III. SETUP FOR DYNAMIC COOLING 
A. System and reservoir 

A reasonably general set-up for dynamic cooling re- 
quires to specify the class of Hamiltonians for the target 
system S and the reservoir R, together with their ini- 
tial state. Here, both subsystems will be taken to be 
finite-level quantum systems with the overall Hamilto- 
nian i/s+R- The initial state of S+R is a canonical 
equilibrium one with density matrix 



-PHs^ 



PS+R 



tr[e-/5^s+R]' 



(5) 



and the initial temperature T = 1//3 (fee = !)■ We shall 
additionally assume that in the initial state the coupling 
between S and R is negligible: 



Hs+R = Hs+ Hn, 



PS+R = PS(S PR. 



(6) 



The action of external fields is described by an interaction 
term in the time-dependent Hamiltonian of S+R: 



Hs+Rit) ^Hs + Hr + Vit), 



(7) 



which is switched on at time t — and switched off at 
time t = t: V{t) = for t < and for i > r {cyclic 
motion of the external field). V{t) acts on both S and R; 
if it acts only on S, no cooling is possible [l^. Hs+R{t) 



generates a unitary operator U that takes the initial state 
(O to the final state of S+R: 

Ps+r{t) = C/ps+Rf/^ (8) 

Ps(7-) == trRps+R(r), Pr{t) ^tisps+Rir), (9) 

where we defined also marginal final states of S and R. 

B. Definition of cooling 

How to define cooling of S in the non-equilibrium set- 
ting? Note that the maximal cooling is always well- 
defined, because it means that the final ground state 
probability of S is [|0) is the ground state of Hs] 



Po(t) = (0|/7s(r)|0) = 1 (or T = 0). 



(10) 



Defining a non- maximal cooling via Po{t) is reasonable 
for two-level systems, since a larger po(''') means that the 
energy distribution in the final state is more shifted to- 
wards the ground state ^. If in addition the state of S 
is diagonal in the energy representation, one can equiva- 
lently express the cooling via temperature defined as in 
(O. For a multi-level system S one needs to take care in 
defining the meaning of a possibly non-equilibrium state 
of S being colder than a given equilibrium state ^. 

For simplicity we assume that S is a two-level system 
with energies * 



< e, 
and equilibrium probabilities 

po = (0|ps(0)|0) ^ 



-13c 



Pi 



(11) 



fl2) 



1 + e-^^' "" l + e-0^' 
The reservoir R is a M-level system with energies 

= Mo < Ml < • • • < PM-i ^ M, (13) 



® One may define cooling via the maximal eigenvalue of the den- 
sity matrix, and not the ground-state probability. In the optimal 
regime of our set-up both definitions agree with each other, be- 
cause the final state of S is energy-diagonal; see section FlV Al 

"^ When comparing two systems with the same energy levels — e.g., 
the target of cooling before and after the cooling process realized 
via a cyclically changing Hamiltonian — one can define cooling by 
requiring that the whole energy distribution is shifted towards 
the ground state. In effect, this amounts to using majorization as 
measure of cooling; see 28] and Appendix A for the definition of 
this concept, see also [21] for a related approach. The drawback 
(or viewing differently an advantage) of this definition is that 
when the number of energy levels is larger than two, not every 
two states can be said to be cold or hot relative to each other. 
Another (less preferred by us) approach to cooling would be to 
employ global measures such as entropy. 

* Consider a multi-level system whose energy gap between the 
ground-state and the first excited state is smaller than the gap 
between the first and second excited state. For sufficiently low 
initial temperatures this system can be regarded as effectively 
two-level system, because populations of the second and higher 
energy levels can be neglected. 



and initial, equilibrium probabilities [see ^ [6])] 



o-l^l-l-l 



TT; 



i + Efei7'e-'""= 



l = 0,...,M-l. (14) 



The eigenvalues of the initial density matrix ps+R read 



{^fc}fc=0 ^ = (P07I"0, PlTTo, P07I"1, PlTTl, . . .). (15) 



C. Unattainability of the absolute zero 

The above set-up suffices for showing the unattainabil- 
ity of the absolute zero, po{t) > 1, for the target of cool- 
ing S, given the initial state ([5]) of the reservoir. The ffiial 
ground-state probability po (r) of S follows from ([5l- [T?|l 



,M-1 



. nE P.7r^[l-|(0,a|C/|z,7)P],(16) 

where {|j)}i=o ^^^^^ {|Q^)}fio^ ^-''^ the eigenbases of, re- 
spectively, Hs and Hji in ([6]), and {tTq}*^^ is given 
by (HI. Since ELoE*lro'l(0,«|f/N, 7)1' = 1, some in- 
equlities 1 > |(0, a|C/|i,7)p must be strict. Noting also 
PiTT^ > we deduce from (|16p 



Po(''") < 1 (strict inequality). 



(17) 



This argument applies for M = cxd, where, e.g., 

Ei=oE7=ql(0'"l^l«'7)P = 1 and E^=o^7 = 1 are con- 
vergent series, and where U is unitary for M = oo. Note 
that the argument leading to ((T7| is essentially based on 
PiTTj > and hence n^ > 0; see (fT4|) . 

With trivial changes the derivation of ([TTl) applies for 
a multi- level system S. 



D. Lower bounds for work 

Here we summarize restrictions imposed by the second 

law on the work needed for cooling^ ^Our presentation 

uses the same ideas and methods as [33, Isj] ; see also [3a| 
in this context. 

Recall that for the considered unitary (thermally iso- 
lated) process the work done on the system is equal to 
the [average] energy difference 33, 34]: 



W = tr( [ps+r(t) - Ps+ii]Hs+-r) 



(18) 



where we recall that the interaction V{t) with external 
fields is switched off after the final time; see ([7]) . 

Since the initial state ([5]) is at a Gibbsian equilibrium, 
any unitary operator that changes this state costs some 
work. Indeed, given the initial state ([5]) and the unitary 
dynamics implemented by a cyclically changing Hamil- 
tonian ([71 [8]), the work ([T8|) invested in S+R amounts 
to 



pw 



= tr(ps+R In ps+R - Ps+r(t) hi ps+n) (19) 

= tr(ps+R(T) In Ps-i-r(t) - Ps+r(t) In ps+r)(20) 
= 5'[ps+R(T)llps-t-R] >0, 



(21) 



where in moving from ([T9)) to ([20| we used the uni- 
tarity of dynamics; see ([8]). The relative entropy 
'S'[ps+r(''')11ps+r] is non-ne gati ve and nullifies only if 
Ps-i-r(''') — Ps+r; see Refs. [53] for further features of 
the relative entropy and its role in (quantum) informa- 
tion theory. Thus any change of the initially equilibrium 
state via a cyclic Hamiltonian process costs some work. 
This work is also a resource of dynamic cooling. 

In the context of cooling one can derive more stringent 
lower bounds on the work. Using ([6]) we rearrange ([20]) : 



PW = ^[ps(r)||ps] + 5[pr(t)||pr] + /sR(r), (22) 

^s+r(t) = tr[ps+R(T)lnps+R(T)] 
-tr[ PS (t) In PS (r) ] - tr[ pR(r) In pR(r) ] > 0. (23) 

/s+r(t) is the mutual information between S and R in 
the final state. /s+r(''') > due to the sub-additivity of 
the entropy. Now ([ll IH [22|) and ([HI [23|) produce 



W > r5[ps(T)||ps] 

= F[ps(r)]-F[ps] = A^>0, 



(24) 



where for any density matrix a, the free energy at the 
initial temperature T is defined as 



F[a] = tr(iJscr) + Ttr(crlncr). 



(25) 



If cr is a Gibbsian density matrix at temperature T, F[a] 
coincides with equilibrium free energy [cf. ([3])]. The 
meaning of ([22]) is that the work equals the sum of the 
free energies differences for S and R + the energy stored 
in the mutual information between S and R. 

AF in ([24)1 is different from the free energy difference 
([4]) of the Nernst set-up, where both the initial and final 
states are at equilibrium. Hence in ([4]) the temperatures 
Tin and Tgn refer to the initial and final states, respec- 
tively. For dynamic cooling only the initial states are at 
Gibbsian equilibrium ([5]). Hence the free energy in ([24l 
I25p refers to the initial temperature only. To compare 
quantitatively ([24)1 with ([4|) , we note that g-m in ([JJ refers 
within dynamic cooling to the S— R (system-reservoir) in- 
teraction that is zero both initially and finally. Assume 
that the final state ps (t) of S is Gibbsian (canonic equi- 
librium) at temperature Tfin. Then the thermodynamic 
entropy in ([3]) is the von Neumann entropy from (j25l) . 
Recalling that Tfin < Ti„ = T we obtain from ([H [24|) 

F[Tan,gin] - F[Tin,gin]-TS[ps{T)\\ps] 

= -tr[ps(r)lnps(r)] (Tfin - Tin) > 0. (26) 

Thus the difference between ([M)) and ^ tends to disap- 
pear with the final entropy of S. Such a situation will be 
met in section [VH 



E. Summary of conditions 



In analogy to section lll Dl we summarize the dynamical 
cooling set-up in terms of the following task list: 



Dl. Initial state: canonical equilibrium system and 
reservoir; they are uncoupled. 

D2. System (target of cooling) : microscopic. 

D3. Cooling dynamics: unitary in the space of the 
system + reservoir. Not necessarily quasi-equilibrium. 

D4. Reservoir: possibly microscopic. Not necessarily 
larger than the system. 

D5. System-reservoir interaction (driven by external 
field): changes cyclically. 

D2 makes obsolete the need for the entropic formula- 
tion of the third law, since a thermally isolated process 
done on a finite S is not uniquely characterized by its 
entropy [32] (instead the full spectrum of the density ma- 
trix determines the set of states that can be reached by 
unitary processes). In this paper we assumed D2* rather 
than D2, i.e. that the system S is two-dimensional. This 
more restrictive condition should be relaxed in future 
studies. D5 is motivated in the same way as N5. To 
grasp the difference between N4 and D4, we shall study 
in section rvTl the dynamic cooling set-up under condition 
N4. 

As for the tractability we keep the conditions Tl— T3 
from section IIIDl The unattainability of the absolute 
zero is proven in section UlI CI via Dl and D3 — without 
explicit reference to Tl— T3. Nevertheless, T2 will be 
needed below for understanding limitations on reaching 
lowest (non-zero) temperatures, while T3 is inherently 
demanded for applications of the dynamical cooling: The 
cycle time now refers to the time needed to realize the 
unitary operator ([7]) via a suitable Hamiltonian. 



IV. OPTIMAL RESERVOIR 



The task list of section IIII El does not yet fully specify 
the reservoir. In this section and the following two we 
will investigate three different models; their respective 
parameters are then shown to characterize the way in 
which cooling can be achieved. 



A. Max-Min cooling scenario 

The purpose of dynamic cooling is to increase the 
ground state probability of S subject to the constraints 
listed in section [III El Obviously, some set-ups will work 
better for the desired task than others. For studying 
principle limitations we should thus optimize the design, 
i.e. try to maximize this very probability. This has to 
be done i) over all unitary transformations C/ in ([8]); ii) 
over the energy level distribution of the reservoir assum- 
ing that they are bound from above, i.e., /x in (J13p is a 
finite, fixed number. This parameter of the upper energy 
level of the reservoir will turn out below to control dy- 
namic cooling. We denote these maximization strategies 
as max[/ and max^^. respectively. The overall maximiza- 
tion is max^/^^j. . 



Generally, the above maximization procedure will 
specify a set of equivalent scenarios only, i.e. leave some 
free parameters (see below). Additional requirements 
may be imposed to reduce this ambiguity. Here we at- 
tempt to minimize the work needed for given cooling ef- 
fect. This is why our intended goal requires a "max-min" 
scenario: first the ground state probability of S is maxi- 
mized and only after that the work is minimized. 



B. Maximal cooling 



Assume that R has an even number of energy levels: 



M = 2n. 



(27) 



Appendix [X] shows that the unitary operator, which 
leads to the largest final ground-state probability for S, 
amounts to permuting the elements (jlSp of the initial 
density matrix ps+n- We recall that it is diagonal in 
the energy representation. Hence in the maximal cooling 
regime the final state of the two-level target is diagonal 
in the energy representation: it has a well-defined tem- 
perature, which is hence not imposed, but emerges out 
of optimization. 

The fact that the two-level target of cooling ends up in 
an energy diagonal density matrix (with the ground-level 
probability greater than the excited level probability) is 
straightforward to establish: otherwise, there will be a 
unitary operator acting only on the two-level system such 
that its ground-state probability is increased. 

Consider the eigenvalue vector of the final state of 
S+R with the largest ground-state probability po{t) of 
S. In this vector the largest 2n elements of the vector 
([T5|) are at the odd places [counting starts with 1], and 
the final, maximized ground-state probability for S reads 



ma.xu[po{T)] = smax2„[{wfc}^"o^] 



(28) 



where smax/j [a] returns the sum of k largest elements of 
vector a. Finding the 2n maximal elements of {wfe}j,"g^ 
leaves some freedom in U. One represents (l28l) as 



En— 1 
TTk+PlTTo 
k—0 

-fsmax„_i[pi7ri, . . . ,pi7r„_i,po7r„, . . . ,PoTT2n-2], (29) 

where we used (jlSp and straightforward induction over 
n. 

Maximizing (|^5|) over the energy levels (fT5)) amounts 
to maximizing it over the Boltzmann weights Vk — 6"^^*= 
[see p^fH]) ] under constraints 



1 > «i > ■■• > ^^2n-l >v = e 



'/3m 



(30) 



Since the LHS of (|28p is a ratio of linear functions of Vk, 
it can maximize only at the borders of the allowed range 
of Vk, i.e., some of w^'s are equal to 1, while others are 
equal to V. The physical meaning of this result is that the 
optimal reservoir is an effective two-level system, a fact 



that greatly simplifies searching for the maximal (final) 
ground-state probability of S. To illustrate this result 
consider an example in ([29l) : n — 2, sniaxi[pi7ri,po7''2] = 
max[pi7ri,po7''2] — PiT^i- This reduces ([2^ to ttq + tti, 
which maximizes for ui = 1,^2 = w. 

It now remains to check all possible arrangements of 
energy levels that render R an effectively two-level sys- 
tem. This check produces the maximum of (j29p for 



fJ-o 



fJ-n-l = 0, 



fJ'T, 



^J■2n-l = ^J■■ (31) 



Hence the optimal reservoir for cooling a two-level system 
has to have a degenerate ground-state for n > 2. The 
maximal ground-state probability now reads 



max;7,^Jpo(T)] = 



E 



1 



VTfe = — 
fc=0 1 



-/3m 



(32) 



-e/T„, 



where 



Tlnin = Te/fi. 



(33) 



Thus (j32p is the maximal probability for a reservoir with 
maximal energy /i; Tmin is the minimal temperature. 
Note that ([55]) is consistent with the unattainability ar- 
gument (|17p. For /Lt — >■ oo, Tinin would approach zero. 

We get cooling, i.e. max(7^^^[po(''')] > Po, only when 
the largest energy of R is larger than the energy of S: 



fj, > e. 



(34) 



This asymmetry ensures that the two-level is cooled in 
presence of the reservoir. 



C. Minimal work 

We turn to calculating the minimal work given the 
maximal probability ([5^. Write down p3|) under ([5T|): 



P0TT,pin, . . . ,poTT,pin,poTT,piTT, . . . , PqTT, piTT, (35) 



2n elements 



2n elements 



and TT 



where we defined tt = — tt-i — r o,iiu /. — —n-: — r- 

n{l+v} n[l+v) 

Eq. ([5^ is found after the first 2n elements in psp are 
distributed over the odd places in the eigenvalue list of 
pI+r [counting starts from 1]. Concrete places they oc- 
cupy are not important for obtaining the optimal cooling 
effect (j32p : Those various orderings constitute an equiv- 
alence class. Concrete places, however, become impor- 
tant for minimizing the work. Note that for a fixed final 
ground-state probability of the two- level system minimiz- 
ing the work amounts to maximizing the overall proba- 
bility of the lowest level ^o — . . . — fin-i — for the final 
state of the reservoir; see (I5T]) . Then (|55|) implies that the 
minimal work is obtained for a permutation that does not 
touch the underlined elements in (1351) . but interchanges 



those that are not underlined: each piir with some el- 
ement pitt; see Appendix |^ for the argument reducing 
the considered unitary to a permutation. Recalling that 
the initial energy of the system4- reservoir is ff _^^ +Pi£ 
and using ([T5)) we get 



W^{f,-e) 



1 



-/3m 



PO 



(36) 



for the work. Now ((36|) shows that reaching the lowest 
possible temperature T^-^ requires the work 



W C^i UPl r^ 1/Trnin- 



(37) 



This constitutes a parametrization of the attainability 
constraints for dynamical cooling: Work W and and the 
inverse temperature reached, l/Tmin, are proportional. 
Accepting that W cannot be infinite (see condition T2 
in section IIIDI and IIIIEp , we have to accept that the 
minimal temperature reached cannot be zero. 

Let now the number of energy levels of R be an odd 
number: AI = 2n+ I. Instead of ([29]) we get 

En— 1 

+smax„[pi7ri, . . . ,pi7r„_i,po7r„, . . . ,Po'iT2n-i]- (38) 
The final ground state probability is maximized for 

fiO = ... = fln-l =0, ^J.n = ■■■ = ^J■2n = Mi (39) 

meaning again that the optimal reservoir has to have a 
degenerate ground-state for n>2. The maximal ground 
state probability reads: 



niax;7^^j^ 



n+p^e ^^^ 



n+poe 



■-/Sm' 



(40) 



Condition (|34)) is still needed, and the qualitative conclu- 
sion from studying the work cost is the same as above. 

In constrast to ([5^. expression (HUl) already depends 
on the initial probability po of S. Eqs. (|32l) is recovered 
from (|40|) for a many- level reservoir n ^ 1. Instead of as- 
suming such a many level reservoir with the optimal uni- 
tary operating on the joint Hilbert space of this reservoir 
and the target system, we can relate (|40)) to (l32t under 
a weaker condition. Apply the cooling protocol repeat- 
edly with the reservoir — having a finite, odd number of 
energy levels — re- prepared in its equilibrium state, e.g., 
via fast relaxation, as it happens with electronic spins in 
solid state NMR Q. Then the ground-state probability 
increases iteratively as [see (j40|)]: 



['-HI] 



+ p[;ie-^'' 



1 = 1,2,. 



' n+{n+ l)e-l^>'' 
For / ^ 1 the result of this iteration converges to ([5^ 



(41) 
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D. Trade-ofF between maximizing cooling and 
minimizing work 

The existence of the above work cost for dynamic cool- 
ing raises the following question: to what extent can we 
reduce this cost if, given the upper bound on the reser- 
voir energy spectrum [see (|13|) ]. we sacrifice some ground- 
state probability of the target system, i.e., instead of 
reaching ([5^ we agree to reach a somewhat lower final 
probability. As compared to (|32]) . this will require a dif- 
ferent reservoir and a different unitary transformation for 
cooling. The answer to this question is that even a small 
decrease in the achieved ground-state probability can sig- 
nificantly reduce the work-cost. We shall illustrate this 
fact via an example. 

Let us take M = 2n = 4, and postulate a unitary 
operator that permutes the eigenvalues of the initial state 
(fT5|) so that the final-state eigenvalues read 

(PQTTO, P07r2, PlTTo, PlTT2, PqTTi, PqTTs, PlTTl, PlTTs). (42) 

The logics of obtaining (|42l) from P3|) is as follows. The 
first four probabilities pqttq, pittq, Pqt^i, Pit^i in (|15l) are 
moved to odd places producing 

Po(r)^^o + 7ri= ^ , _,.^^ ^ _„^^ , _., ■ (43) 



1 



-fit^i 



~Pf^2 



-fSt^s 



which is similar to ttq -I- tti in p2p . The remaining four 
probabilities paTT2, Pitt2, Pot^s, Pit^s [last four elements 
in (|15p ] are then arranged in between without changing 
their mutual order. The work and the final probabilities 
for the four reservoir energy levels read 

W = e{pa-poiT)) + 51'Li'"'^'"'^^^ ~ '"')' (^^^ 
TToir) =_po(7i'o + 7i'2), 7i'i(r) =pi(7ro + 7r2),(45) 

7r2(T) = pai-Ki + TTa), 7r3(r) = pii^i + 7r3).(46) 

Let us assume that 



/3m2 =(3fi3 = f3fi:>l 



(47) 



is a fixed number. Hence /ii is the only parameter over 
which we can minimize (j44p . In the limit /3/i — >■ oo both 
and P5)) converge to 1: 



max,7,^Jpo(T")] 



1, 



Po(t-) 



1, 



while their ratio is a finite number 



l-maxc/,^Jpo(T)] _ 1 „^ ^^ 



(48) 



(49) 



where C'[e~''^] can be neglected due to (H71) . 

We now show that although (|32|) and (|43|l are quite 
close to each other, the work needed to obtain (P5| devi- 
ates significantly from ([36)) . Note that 7r2 and tts are ex- 
ponentially small with /3/x — >■ cx). We neglect such terms 
in (l44l) and obtain 



M^ 



^epi -|-^i(pi7ro -TTi) +/i7ri, (50) 



/3M^ = pi In 



Pi 



/3mi(pi 



-/3a*i 



) -h /3/ie-'3^i 



1 



Pi 



-/9a«i 



(51) 



To get dnH [so]) from (|44H46l) we employed p2l IM)) . It 
should now be clear that for /3/i ^ 1 we get a logarith- 
mically growing work, PW = Pi lii[ j'i ], if we choose 
/3/ii = ln(/?/i). The minimization of the RHS of (|5T|) over 
/3/ii satisfying /ii < /i produces a very similar result. 
For /3/i ^ 1 neglect in (I5T|) all terms with e"^'^^ except 
the factor (3iie~^'^^, where a small term e^'^^i is multi- 
plied by the large f3fi. Now substitute /3/ii 
into (j5ip . where ?/ is a new variable 



PW = piln 



+ pi{y + In 



,3m 



In i^ 

IVPl 

This produces 
After differen- 



tiating over y we obtain for the minimum; y = 1, or 



PW = pi In 



1-pi 



(52) 



and thus a logarithmic rather than linear increase of W 
with /i, i.e. with /3'. In the limit /3/i — >■ oo the correspond- 
ing ground state probability is only slightly reduced: 



l-maxc/,pJpo(T)] _ 1 
1-Po(t) "2 



1 



P]_ 

/3^ 



(53) 



This equation can be rewritten in terms of tempera- 
tures. Recall (|481 l32t . and introduce a temperature T 
via po{t) = (1 + e-'/^)-i. Then dSS]) reads 



1 



1 



T T„ 






2 V /3m 



(54) 



For /3/i — > oo the r.h.s. of this equation goes to — 1£^. 

If one allows further deviations from the maximal prob- 
ability (|32|) [more than (|53|) ] , then it is possible to reduce 
the work even below the logarithmic dependence (|52p . 
Appendix |B] shows that the example (|42l) generalized 
to n = 4 (eight-level reservoir) produces [for /3^ ^ 1] 
the double-logarithmic scaling of the minimimal work, 
W ~ ln[ln[/3/i]] provided that the deviation from the 
maximal probability scales as 



1 -max[/,^^[po(r)] 
l-Po(r) 



i<> 



O 



1 



HPfi) 



)• (55) 



Realistic reservoirs tend to deviate from the optimal 
model studied above. In the following we investigate 
models subject to additional constraints. 



V. HOMOGENEOUS SPECTRUM OF 
RESERVOIR 

A. Finite number of levels 

We saw above that the maximal energy ^ [see (|13p ] is 
an essential parameter for dynamical cooling via a reser- 
voir with a finite number of energy levels: cooling to 
ground-state probability 1 is possible only for ^ — >■ oo. In 
the optimal scenario p2p (with two levels only) fi conin- 
cides with the energy gap. 



Here we study the simplest cooling scenario, where 
/i can be large without increasing the level spacing. 
This scenario also illustrates the limit of the inifnite- 
dimensional Hlibert space for the reservoir. 

We assume that R has M equidistant energy levels: 



M/c 



Sk, 



0,1,. 



,M-1, 



(56) 



where <5 > is the reservoir energy gap. We realize 
with respect to this reservoir the same max-min scenario. 
However, no optimization over the energy gaps of the 
reservoir will be carried out, i.e., 6 is a fixed parameter. 
The initial-state eigenvalues of S+R read: 






,p,v''~\p,v''^-^\ 



(57) 



,M-1 



where we defined 

V = e-P\ Z = EH7«' = (1 - ^'')/(l - «)• (58) 

We define a from 

ln(po/pi) 



Pov" >Pi> Pov" 



ln(l/u) 



(59) 



where \_x\ is the floor integer part of x, e.g., [0.99J = 
0. Let us re-order the elements of vector ([57)1 (first the 
largest element, then next to the largest, etc): 



<^ = -^[ Po.Pqv, 



,Pov 



Pi,Pov°'~^^,piv,pav°''^'^, 



Q + l 

M-Q-2 A/-1 

,piV ,poV 



2(M-a-l) 

pi«*^— ^pX^-^...,Pl^;^^-l■ 
^■ V ' 

a+1 



(60) 
(61) 
(62) 



where the curly bracket shows the number of elements in 
each group. In (|60l - f62|) we select the first M elements; 
their sum will give the maximal final ground-state prob- 
ability of the target two-level system. Then the optimal 
unitary transformation amounts to distributing those M 
largest elements over the odd places in the eigenvalue list 
of the final density matrix. Clearly, M = a + 1 means 
that no cooling is possible with the considered reservoir. 
Let M — (a + 1) be an even number: 



M- (a + 1) = 2s, 
where s is an integer. Now (|61l) is factorized as 



(63) 



Thus the final ground-state probability of the target two- 
level system reads ^ 



1 



Q + S 



maxc/[po(T)] = ^[PoJ]^^^!''^ +Pi^^^^«'=] 

1 — piw* — Pqv"^''^^ 
= 1 - vM ■ 



(65) 



As expected, (|65|) is smaller than the bound (|32p with 
fi = S{M- 1) . 



B. Infinite number of levels 

For finite S and a we get from (|65p : 

ma,xu[po{T)] ^>- I when i\/ — > oo and s — > oo. (66) 

The minimal work necessary for (j66l) is finite. In show- 
ing this let us restrict ourselves with the simplest case 
a = in (gS]). Eqs. ((60ti62)) now read 



po,p^,...,p^v(''-'^/^pov^^'-'y^ 



piv^''~'y',Pov^''+'^/^ 



M 

,...,Pov''-\p^v''~' (67) 



M 



Eq. (j67p implies for the minimal work in the limit (I66II 

^ ^^ ^ (M-3)/2 

W^e{po-po{T)) + -ipi+Pov^^) Yl (2^ + !)«'' 

fc=0 
(M-3)/2 M-1 

+ ^{pov+piv^^) Yl (2/c + 2)v^ - - ^ fcv'= 



fc=0 



k=0 



= {S- e)pi + 



Sv 



(68) 



which is clearly finite, positive and is larger than the 
bound AF = -i ln(l + e"^'') demanded by ((24)) . 

The convergence ([M)) does not mean that for 
Af = oo (harmonic oscillator spectrum) we would get 
rB.axij[pQ{T)] — 1 (or Tmin — 0), because for A/ = oo the 
transformation that leads to (|65p is not even bijective, 
let alone unitary (i.e., the limit of unitary processes for 
Af — >■ oo is not unitary); see in this context our discus- 
sion after (IT7|) . Let us take a = in (|5DH5^ . For a finite 
M the optimal unitary — which we recall amounts to a 
permutation — distributes the first half of the vector ([60]- 
[62| over the odd places in the final vector. The elements 
from the second half are distributed into even places of 
the final vector. For Af = oo this second (sub)process 



6Tj) =Pi,PqV 



s — 1 a-\-s 

.,piv ,pov ^ , 



Piv ,pov 



a+s+l 



Af-Q-1 



M-a-1 



® Likewise, for an odd M ~(a + l) we introduce Al — a — l = 2u + l, 
(64) where u is an integer. Instead of II65I I we get ina,xi/[po(T)] = 
~P^^ ~Po^ ^ which implies the same consequence I I66II . 



1_„M 
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disappears in infinity. Hence the limit of tlie above per- 
mutation for M — )> cxD is not even bijective. 

Cooling with the harmonic oscillator reservoir (Af = 
oo) is a well-defined problem provided that one ensures 
that the operator U stays unitary for AI — oo. For each 
such unitary we have Poi^) < 1; see (|17|) . But the maxi- 
mum maxu[po{T)] now does not exist. It can be substi- 
tuted by supremum 



suPc/[po(t)] = 1, 



(69) 



as follows from (|66p . Eq. (IMl) points out on an important 
difference between the infinite and finite-level situation. 



For simplicity we focus on the even a case ((7T|) . Then 
taking in (|BDH5^ all the initial elements (together with 
their degeneracies) up to the index m means dividing 
it into two equal parts. Thus the largest ground-state 
probability achievable is equal to the sum of the elements 
from the first half of uj: 



max[/[po(T)] 



1 



[l + v] 



N 









a 
2' 



Examples of (17^ are presented in Fig. [51 



(72) 



VI. RESERVOIR CONSISTING OF N 
IDENTICAL SPINS. 

We study this modular case for three reasons: First, 
in the thermodynamic limit A^ — > oo the reservoir will 
become a standard thermal bath; hence one expects to 
establish connections between the Nernst set-up [section 
HI] and the dynamic cooling set-up [section IIIIj . Sec- 
ond, the model is relevant for polarization transfer exper- 
iments, and thus was widely studied — albeit for the high- 
temperature limit only — in the NMR literature [111, UM ■ 




P 



A. Maximal cooling 

Now reservoir R consists of N identical spins, each 
one with energies and 5 > 0. The energy levels of the 
reservoir, Ea = 5A^ are equidistant with degeneracy 



dA 



N\ 



[N ~A)\A\ 



(70) 



One can apply ([SOHS^ to this situation — with a being 
defined as in ([5^ — but now each element piV^ has to be 
repeated dA times. The indices under the curly brackets 
in (j60H62p indicate now the number of distinct elements. 
Again, we need to divide the vector w in (ISUHS^ into 
two equal parts. For simplicity [and without altering the 
asymptotic formulas (|751 [76)) ] assume that N is an odd 
number. 

If a is an even number, dividing u into two equal parts 
amounts to finding an integer m that satisfies 






k= 



dn 



^--^k=0 



d,. 



,JV 



(71) 



since now the first half of cJ ends up with element piu™. 
This leads to m = [-^J — ^, where [a;J is defined after 

If a is an odd number, the first half of Co ends up with 
element p^v"'^'"'^, and m = [yj — [^J is found from 



E 






fe=0 



dfc = 2 



N 



FIG. 2; The relative cooling ■"'^''^[po(r)]-Po calculated accord- 

° PO 

ing to (|12ll72p versus the initial inverse temperature /3 = 1/T 
for various values of e and A'^ = 201, 5 = 0.1. The relative 
cooling is maximal at some intermediate /9(e). The maximum 
is sharper for smaller values of e. 

Let us now take the limit iV ^ 1. Despite this limit 
both e and 5 are finite; we neglect a in ([72]) by taking 
there m = \J^\ [recall that a = for e < (5, an almost 
degenerate case] and transform ([7^ via 



(l-hw)^^^fc=o 



dkv^ = 



Efe=o dkv^ - E 



, , 1 dkV 

fc— rn+l ^ 



(1 + z;) 



N 



m+lY^^"™"! fc rffc+m+1 



[l+v]'^ Z^fc^O 



(73) 



Since m = [-jj > Nv/{l + v), v'^dk+m^i is a decaying 

is dominated by 



2 

function of fc; hence X]fc^™~ ^'^ k+m+i 

its first few terms and is bounded from above for A^ — > oo 



We are left with estimating — 



' + 'd„ 



(1+^ 



for TO 



N 



For numerics one can use the hyper-geometric function: 



E 



JV-m-l fc gfc + m+l 



1+1 



2-Fi [1, m - A^ + 1, 2 -I- m, -?;] that holds 



for any N — m — 1 > 0. This formula is derived by express- 
'^iil. via Euler's gamma functions [recall I I70I1 1 , extending 

= Efclo' a°<i employ- 



ing 

the summation to infinity: Ylk=o^ 

ing the standard definition 2F\{a, b, c; z) = J^^ 



=0 (c)fc 



fc! 
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Using the Stirling's formula (|74|) 

N\ ~ ^/27rN{N/e)^, 
we get for A^ ^ 1 



where 



(74) 



lln(l-maxt;bo(r)]) = -5[i/||a]+0(^), (75) 

where the relative entropy Sf^/Hcr] between the proba- 
bility vectors (5,5) and (y^, j^) = cr satisfies 

5[-7||a] = In ^-^ = Incoshf^l; (76) 

see also (HJ). Note that a is the state of the single reser- 
voir spin [the reservoir consists of N such identical spins] . 
Note that for A^ ^ 1 the asymptotic behaviour of po (''') 
does not depend on e. 

Thus for A'^ — >■ 00, the (small) deviation of po(''") from 
1 is controlled by the relative entropy (1751) . The final 
temperature of the two-level system reads from ()751 I76p 



^ min — 



NS{\l\\a\ 



(77) 



Note that (j77p is consistent with the unattainability ar- 
gument (fT7| . 

Let us now turn to calculating the minimal work 
needed for this cooling. 



B. The minimal w^ork 

The message of the following calculations is to show 
that the minimal work required for cooling (|751 1761) can 
converge to its thermodynamic lower bound (j24l) in the 
thermodynamic limit A^ ^ 1. 

We start by the case e = 0. Conceptually, this case 
of degenerate target system energy levels is interesting, 
because no brute force method (i.e., a low-temperature 
bath) can ever increase the occupation of one of the en- 
ergy levels. In addition, the work done for cooling is not 
blurred by the energy released from the two-level system. 
Starting with this case is also useful for technical reasons, 
moreover that e > will not bring essential news. 

As follows from the discussion around (|71l) , the vector 
can be written via its two halves as 



(a, 6), 



(78) 



of the hyper-geometric function 1591 . Here {a)j. = y{\ ^^'^ 
(— a)fc = (—1)* p, _fc I 1^ are the Pochhaniiner symbols, and r(2) 
denotes the gamma function. 



1 



2(1 +t))^ 

1 

2(1 + t;)^ 

c=[0,l,2,...,7V], 



z;0,i,o,...,«[^/2l,t;[^/2l], (79) 



(81) 



and where c is the vector of the reservoir energies. Recall 
that in (175)) and ([50]) each factor v^ is repeated d^ times. 
Likewise, in ()8ip each symbol k is repeated dk times. 
Vectors a, h and c have equal number of components 2''^. 
The discussion around (1711) implies that a contains the 
largest 2^ elements of the vector ([60l - [62|) . Thus the min- 
imal final reservoir energy needed for the optimal cooling 
([75]) is the following inner product 



8{d+h)c. 



(82) 



Once reservoir's initial energy 5Nv/{l + v) is known, the 
initial and final energies of the target spin are known as 
well, the minimal work is determined via (|82p. 
Note that dkv'' is peaked around 



A^- 



NL 



(83) 



where \x\ is defined after ([55]) . v = e ^^ and the dif- 



ference 



•» l+v 



fluctuates for different A^ 3> 1 only by 

the amount 1/A^. Hence in a we can focus on dominant 
energy levels 



{<5(A^^ 



= -A' 



(84) 



where A is determined by requiring that the ratio 
"rf^^t)"g — '^^ ^^^ central (maximal) probability to the 
probability at the edge of the interval ([M)) is exponen- 
tially (over N) small. This guarantees that considering 
only (|M)) suffices for calculating quantities that are not 
exponentially small. Employing for A^ ^ 1 Stirling's for- 
mula ([74|) . we get 



JiV^iA 



,,7VC±A 



dNe, V 



Ni 



— fi 2£(1-£)JV 



'(^)+o(i 



■) 



(85) 



For achieving the sought exponential smallness we need 
to require A = N°'~^^''^, where a is a fixed small number, 
e.g., a = 1/10. Hence the total probability of energy 
levels that do not fafl in the interval dHU) is 0[Ne^^^''] -^ 
0. Now in (j82|) we can neglect 6c and keep in ac only those 
elements of a whose energies lay in the interval (|M)) . 

Within interval ((84|) we select a segment consisting 
of ti™+^5's only. In a it occupies positions with num- 
bers from 1 + J2k=o"^~'^ 2dfc to X^fcio"™ ^dfc- Its length 
is 2dm+Ni- It appears that there are only two energies 
^m — 1 and $m that correspond to that segment in c. 



We write $m = NS, + m 
integer that satisfies 



£m, where im is the minimal 



dk 

k " 



'k=0 



2dk = D„ 

k=0 



> 0. (86) 



The reason of having only two energies for each segment 
in the dominant interval (|M| is that £„ does not depend 
on m, as we see below. 

Each sum in (j86p is dominated by its largest summand. 
Recalling from ([7D|) 



d. 



fc+i 
dk 



F^[l + 0(^)]> 



(87) 



N 



we get a general pattern for approximating such sums ^^ 
Eq. ([86)) re-writes as [neglecting factors 0{^)] 






'fc=0 



1-v ' 



(89) 



These relations imply 



/ — / — 



"ln2" 
Ini 

V 


= 


"ln2" 

/3<5 



(90) 



where [x] is the ceiling integer part of x, e.g. [0.99] = 1, 
[1.1] = 2. As anticipated, £m does not depend on m ^^. 
Hence we write for the minimal work (1821) 



Sac = 



'Z 



2{1 + v)'^ '^m=-A 

+ (2djvc+m - A„)(A^e + m + £)] 



where Dm is defined in 



t-^«+"[An(iVC + ^ + ^-l) 

(91) 
Eq. dHU) implies: 



i:*,, 



dm+rn \_^ ' [l+0{-)]. (92) 



Note that 2dN^+m ~ Dm > for the considered range of 



m 
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Using dnH) we get for dHI 
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FIG. 3: The minimal work W{5) for e = and T = 1. The 
dotted line is Tin 2, the lower bound ^^ for W{S). 



TABLE I: The minimal work W as a function of A^ for 5 = 
ln2 = 0.693147, e = and T = 1. 



N 


W 


A- = 601 


0.716891 


A^ = 1001 


0.711592 


N = 1401 


0.708766 


N = 2001 


0.70623 



Neglecting ^^ J^r 



-A (!+,;)« 



we get for the work 



W = d 



vjv-'- - 2) 
2(1 -«) 



(93) 



This expression implies several important conclusions: 

• The initial and final reservoir energy differ from each 
other by a factor W = 0(1); hence the state of each 
reservoir particle changes by a quantity of order 0{1/N), 
neglegible for TV :^ 1. 

• When -§T- approaches to an integer number from 



fS5 



below, dnniinni) predict: 



W^^Tln2. 



(94) 



According to ((21)) [and recalling that e = 0] Tin 2 is the 
minimal possible work necessary to cool from the abso- 
lutely disordered state [the equilibrium state with e = 0] 
to an almost ordered state (|75|) . Hence the thermody- 
namic bound ([24]) is reachable for certain finite (5's. The 
behavior of W{S) is shown in Fig. [3] For a finite N, W 



^1 Note that both 1(83 and l|88ll contain error 0{^). The error 
0(jj) from H87|l does not accumulate in l|88ll . because the sum 
^l^:2ndk is dominated by few (smaller than A^) terms around 

This fact together with l|90| l and the reasoning above it allows 
to guess an upper bound S£ for the work W. This bound is 
confirmed by I I93I I and allows to deduce quickly the fact of I I94II . 

13 In view of ^2[ this reduces to - > ", ~^ . 

' "^ V 1 — v 



Recall l|70|l and denote dA 



(N-iy. 

(N-A-l)l.A\ 



E/^ md 



— ^^ and using kdi. 



(!+«)' 



Put m = N^+m-N^ 
Afdj._i write it as 



mi: 



A-C 
i=-A-f 



JV^+m 



:,Ni+m 



(1-H^)^ 



EA djv^ + m'" 

m=-A (1 + v)'^ 



N^+m 



where N = N — 1 and § = 1 — ^. Each sum in the curly brackets 
is equal to 1 minus exponentially small terms; cf with I I85II . 
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is larger than Tin 2 even if i#^ is an integer; see Table U 
for numerical results. 

• For 6 > T\n2, ^ yields W{6) = 2(13^ ■ now the 
minimal work monotonically increases with S (reservoir 
energy gap). RecaU from ([75l \76\i that for a large but 
finite TV the deviation oiina.xu[po{T) ] from 1 is controlled 
by 5'[^/||cr], which also grows with 5. Thus enhanced 
cooling demands more work. Moreover, one observes: 



(95) 



WKiAF + TS[p\\a], for T In 2 < (5, 



where AF = Tin 2 is the free-energy difference from ([24| , 
p = ^I [for e = 0] is the initial state of the target two- 
level system, and a is the initial state of any reservoir 
spin. Eq. ([55)1 becomes exact for S ^ Tln2. 

The case e > does not require any new idea as com- 
pared to the above derivation; only algebraic steps are 
more tedious. Hence we quote only the final result for 
the minimal work: 



W 




f — ua 

1+u 



-ff-; for / < 1 

l-\-u ' •' — 



(96) 



l-\-uv "(/ — 1)— liO 
1+u 



-if-; for / > 1 
where v = e^^^ , u — e^""^, a is defined in (j59l) . and where 



ln(l-|-u") 



Ini 



, / = 



v{v ^ — 1 — v") 
1-v ' 



(97) 



For e = both branches of ([96|) are equal, so no condi- 
tioning is needed. 

Our conclusion on the reachability of the thermody- 
namic lower bound (|24|) remains valid also for e > 0: 
consider — — ^-^ reaching an integer number from be- 
low, and I reaching an integer number from above. Then 
(|96| produces W = Tln(l + e~''"). This is the bound (|24|) 
[recall dZlEni)]. 



C. Microcanonical initial state of the reservoir 

So far we considered the reservoir starting its evolu- 
tion from a canonical equilibrium state ([Sj [6|). Another 
notion of equilibrium is provided by the microcanonical 
density matrix, where all the energy levels in the vicinity 
[E — K, E + k] of a given energy E have equal probabil- 
ity, while all other energy levels have zero occupation. 
While the canonical state describes a system prepared 
via a thermal bath, the microcanonical state refers to a 
closed macroscopic system whose energy is fixed up to a 
small uncertainty 2k [23, [30] • Hence k should be large 
enough for [E ~ k, E + k] to incorporate many energy 
levels and account for unavoidable environmental per- 
turbations, but small compared with E |29l[30|. 

The canonical and microcanonical notions of equilib- 
ria are normally equivalent for macroscopic observables of 
systems containing sufficiently many particles; see [36| for 



a review. The microcanonical state can apply for systems 
with a finite number particles provided that the interac- 
tion Hamiltonian does have certain chaoticity features; 
see [29|, [331 and refs. therein. There are situations — such 
as systems in the vicinity of phase transitions, or systems 
with long-range interactions — where the equivalence be- 
tween canonical and microcanonical state is broken even 
for macroscopic observables; see, e.g., [36, 38, 39] and 
[40] for review. Generally, these states are different re- 
garding some important aspects |38l - l40| . In particular, 
the local stability conditions (i.e., stability iwht respect 
to small perturbations) for the microcanonic state are 
less restrictive than for the canonic state (with the same 
temeprature) [lo] . Hence it happens in several physically 
interesting situations that the canonic state is unstable, 
and the microcanonic state is the only description of the 
equilibrium (H-SO]. 

For the reservoir consisting of A^ ^ 1 identical spins 
the microcanonical state concentrated at the energy SK 
with if ^ 1 is especially easy to define: the cIk degener- 
ate energy levels SK [see ([70]) ] have equal probability ^, 
while all other energies have zero probability. Because 
the energy level SK is already exponentially degenerate 
we restricted ourselves with the minimal width k — > 0. 

Given this initial state of the reservoir it is easy to see 
using the same construction as in (|60H62p [see our discus- 
sion after (|70|) ] that the maximal final ground state prob- 
ability for the two-level interacting with such a reservoir 
is just equal to 1: 



maxuipoiT)] 



1. 



(98) 



The result (|98)) does not require A^ — )■ cx) or k ^- 0. It 
is also not specific to the microcanonical density matrix. 
Any density matrix for an M-level reservoir that is diag- 
onal in the energy representation, and that has at least 
[■^J zero eigenvalues would lead to (|M|) : cf. with our dis- 
cussion in section UlI CI For example, it also applies to 
the so called 6'-ensemble, where all the energy levels be- 
low a certain E are equally populated, while the energies 
above E have zero probability. 

Eq. ((98)l was obtained via the standard definition of the 
microcanonical state, but it contains an essential ideal- 
ization: it is assumed that at least [^J energies have 
exactly zero probability. Consider the emergence of the 
microcanonic state. Classically, any single system is in 
a state with definite energy. Weak and inevitable in- 
teractions with environment smear this energy over the 
interval 2k. This however does not ensure that energies 
outside [E ~ K, E + k] have strictly zero probability. 

Quantum mechanically, we cannot know the state of 
a single system unless we prepared it, e.g., by measure- 
ment. One [the only?] general method of preparing a 
microcanonical state is to do a selective measurement of 
energy and then isolate the system ^^. Generally, after 



The selection will produce (Luders' postulate) a density matrix 
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a selective measurement all energies will be populated, 
because the measurement itself may be noisy due to our 
inability to control the measurement interaction and/or 
the system Hamiltonian. Alternatively, the system may 
be weakly coupled to its environment prior to the energy 
measurement, which then necessarily refers to measuring 
the local energy value. For an open system the latter 
observable is not strictly conserved ^^ . Again, both these 
mechanisms will generally populate energies outside the 
interval [E — k,E + k]. It now suffices to give a small, 
but non-zero probability to more than 1^1 energy levels, 
and max[/[po(''") ] in (j98p will be smaller than one, as we 
saw in section IIIICI 

Hence, the unattainability of T = for the micro- 
canonic state of the reservoir is recovered provided the 
microcanonic state contains such tails. Admittedly, we 
currently lack any general condition for how small these 
tails of the microcanonic state could be made. It is also 
unclear which generic [not ad hoc] conditions one has to 
impose on the coupling to environment or on the mea- 
surement noise to ensure the tails needed for unattain- 
ability. In this context one anticipates an additional 
tractability condition T4 — referring to a finite environ- 
mental coupling or a finite measurement noise — to be 
added to conditions T1-T3 in section HTdI and HiTeI ^^. 

Note the difference: the reachability of low tempera- 
tures for a canonical reservoir with a fixed upper energy 
would require an ever increasing energy gap (related to 
an ever increasing work cost), while for a microcanonic 
reservoir (with a fixed upper energy again) reaching low 
temperatures requires small tails around the central en- 
ergy. 

Another essential difference is that there is no work 
cost associated with the microcanonic reservoir: the min- 
imal work necessary for the maximal cooling need not be 
positive (i.e. no work has to be consumed) ; see Fig. S) 
This is related to the general fact that the microcanoni- 
cal ensemble is not passive [3J|: there exist for it a class 
of unitary operations [generated by a suitable cyclically 
changing Hamiltonians] that leads to work-extraction 
[41] . Importantly, the unitary realizing the maximal cool- 
ing for the microcanonical reservoir can be put into that 



W 



oc Tl{E, K)p[nTl{E, k), where pin is the initial state of the system, 
and where Tl{E, k) is the projector on the Hilbert subspacc with 
energies [E — k, E + k]. If pin was a sufficiently smooth function 
of the Hamiltonian (e.g., a Gibbs state), a microcanonic state re- 
sults. If Pin was an arbitrary state we face an additional problem 
of relaxation towards the microcanonic state. 
If the environmental coupling is not weak, the local energy has to 
be properly defined; see 1291 for a concrete proposal (LEMBAS 
principle). 

Conversely, if the unattainability of the absolute zero is regarded 
to be a law of nature, it can constrain the microscopic processes 
of environmental coupling and selective measurement. Albeit 
this possibility contradicts to the current paradigm of deriving 
the laws of thermodynamics from microscopic theories, it needs 
to be taken seriously. 




FIG. 4: The minimal work W (needed for maximal cool- 
ing) versus the initial energy S K of the microcanonical state, 
where the initial state of the A'^-spin reservoir is the mixture 
of j^uN-KV. siisi'gy eigen-states with the energy SK. We con- 
sidered: N = 7 X 10", 5 = 1, e = 0. 



class, as Fig. |4] shows. 



VII. SUMMARY 

In this paper we have studied dynamical cooling of a 
two-level system (as target) in contact to various types 
of quantum reservoirs prepared in various types of initial 
states. Based on operational ideas these schemes have 
been placed within the context of given resource con- 
straints; the operational meaning of the third law can 
thus be clarified. 



A. Canonical initial state. Complementarity 
betweeen cooling and w^ork 

Unattainability of T = is a direct consequence of 
the dynamic cooling set-up; see section ITlI CI In contrast 
to the unattainability formulation of the third law, it 
does not involve unproven or disguised assumptions. For 
a given reservoir with a fixed upper energy starting its 
evolution from a canonical equilibrium, and for a two- 
level system as the target of cooling, we explicitly predict 
the lowest non-zero temperature (j32p achievable within 
dynamic cooling. The lowest temperature ([32]) would 
approach zero whenever the upper energy of the reservoir 
goes to infinity. 

We have been able to re-interpret this functional de- 
pendence in terms of the work intake: reaching the low- 
est temperature demands work growing as the inverse of 
this temperature; see p6|) . However, this linear relation 
only applies for reaching that nominal lowest tempera- 
ture possible: compromising for a slightly increased fi- 
nal target temperature would allow to change the work 
growth regime from the linear to (multiple) logarithmic; 
see section HVDI 
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Next, we studied concrete models of reservoir. For 
them we found the minimal temperature attainable (and 
the minimal work necessary to reach it) by optimizing 
over cooling dynamics. In section I VII we modeled the 
reservoir via a thermal bath consisting of a large num- 
ber, iV ^ 1, of identical spins, and established relations 
between the resources of cooling and certain relative en- 
tropies. The lowest temperature Tmin reachable with 
such a reservoir scales as the inverse of N and the rela- 
tive entropy 5'[i/||a], where / is the 2x2 unit matrix 
and a is the initial state of a bath spin; see ([77]) . The 
minimal work W needed to reach this temperature has 
basically two different expressions. If the reservoir gap 
is sufficiently larger than temperature, 6 > T, we get 
W — >■ AF + S[p\\a], where AF is the free energy differ- 
ence of the cooled two-level system, and where p is its 
initial (hot) state. Recall that AF is the minimal work 
required by thermodynamics for achieving cooling; see 
([M| and (gl). If (5 < Tln2, W is close to AF = Tln2. 
Hence enhancing cooling demands a larger work, but a 
substantial cooling can already be achieved with the min- 
imal work AF demanded by the second law. We should 
thus stress that approaching T = 0, while using only 
the amount of work AF, is feasible for several limits of 
reservoir parameters. 

Note that the existence of the lowest temperature is 
related to a finite upper energy of the reservoir (either it 
really has this feature, or one considers only that class 
of cooling operations, where only finite reservoir energies 
couple to the target of cooling). If the upper energy of 
the reservoir is infinite, the unattainability of T = is 
still there, but the minimum temperature does not exist. 
Instead, the infimum temperature is zero; see section lVBI 



B. Microcanonical initial state 

The unattainability of T = can be recovered also for 
a reservoir starting its evolution from a microcanonical 
state; see section lVI CI However, the origin of the strictly 
non-zero minimal temperature is different: It relates to 
tails of the energy distribution resulting from a weak in- 
teraction with the environment (and thus our inability to 
avoid such marginal couplings). Put differently, the dif- 
ference between canonical and microcanonical reservoir 
concerns the physics of the T — > limit. Low tempera- 
tures in the canonical situation require a large energy gap 
in the reservoir. For the microcanonical case low tem- 
peratures require a sufficiently large system which can 
be prepared in a state with an almost definite energy. 
There is no work cost involved. 



C. Comparison with other approaches 

Ref. [21| studies a cooling set-up with three compo- 
nents: the system S to be cooled, the resource R, and 
the environment E with Hamiltonians iJs, Hn, and H^, 



respectively. S and E start in, respectively, the Gibbs 
states 7s and 7e with the same fixed initial temper- 
ature Tin. Then, [21| allows only for cooling opera- 
tions that are unitary on the Hilbert space of S + R + E 
and commute with the interaction free joint Hamiltonian 
Hs + Hji + He- Hence, no work cost is related to imple- 
menting this unitary and cooling is only possible if the 
resource deviates from its Gibbs state 7r. at the same 
temperature Tm, i.e., a resource being in its Gibbs state 
is "worthless"; see in this context (|22lB^ . If the re- 
source Rn '■= R^^ consists of N identical copies, the 
lowest temperature that can be achieved for S is asymp- 
totically determined by an analogue of (|77l) . where in- 
stead of ^[i/llcr] one has S'[7r||pr] [21|. Here 7r is the 
Gibbs state and pR the actual state of one resource copy. 
This relation between (|77p and the results of ^] relates 
to the fact that in our present setting the reservoir that 
consists of N identical subsystems can (asymptotically) 
be used for cooling if and only if the subsystems are not 
in their maximally mixed state |/; see section lVll Hence, 
relative entropies other than free energies determine the 
value of a resource (reservoir) both in the present paper 
and in [21|. Note that the relative entropy S'[7r||pr] 
must not to be confused with 5'[/0r||7r], which is essen- 
tially the free energy up to constants. As opposed to free 
energy, S'[7r||pr] diverges if pn gets a zero eigenvalue. 
Thus the absolute zero is reachable only for a resource 
state having a zero eigenvalue, an aspect closely related 
to the microcanonic treatment of section IVICI 

Alternatively, the features of the limit T -> can be 
studied via refrigerators that cyclically operate between 
two thermal baths at temperatures Tc and Th {Tc <Th) 
and cool the (finite) low-temperature bath at the expense 
of consuming work from an external source |23l427l | . If the 
refrigerator works at a finite efficiency, then for Tc — > 
the heat taken per unit of time from the low-temperature 
bath scales as r°+i with a > [H, [11], showing that 
cooling the low-temperature bath is progressively slowed 
down. The optimal behavior a = is reached for the 
refrigerator model studied in |27| . 



D. Open issues 

An important aspect of practical implications is the is- 
sue of time required to complete the cooling transforma- 
tion. In our present operational analysis we have dealt 
with this problem rather formally, assuming that once 
the needed unitary transformation is constructed, there 
is always a (time-dependent) Hamiltonian that realizes 
it. Provided that we do not restrict the magnitude of the 
external fields realizing the Hamiltonian, this realization 
could, in principle, take an arbitrary short time, i.e., the 
external fields can be assumed to function in the pulsed 
regime. However, if we demand that the optimal cool- 
ing unitary transformations are implemented in terms of 
a well-defined base of realizable unitary operations, the 
cooling process may turn out to be complex: it may take 
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a long sequence of the realizable operations to contruct 
the needed unitary. This complexity (and hence time) 
resource needs further studies. 

Another open issue concerns target systems with more 
than two energy levels. This should be interesting es- 
pecially with respect to non-equilibrium aspects. More- 
over, cooling to very low temperatures may in fact de- 
mand considering many-body targets of cooling, because 
at such low temperatures the standard assumptions of the 
weak coupling between the target and the reservoir — as 
well as between different parts of the target — may be bro- 



ken [58|- We nevertheless expect that the basic features 
of dynamic cooling discussed here will survive for more 
general models (the unattainability of the ground state 
and a non-trivial work cost for cooling). 
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Appendix A 

Here we show that the unitaries implementing opti- 
mal cooling can be chosen to be permutations of energy 
eigenstates. 

We are given two sequences A = (Ai, ..., Xn) and D = 

(di,...,d7v) of real numbers. Let {Xl}^^i and {rftl^i 
be the non-increasing arrangements of their elements: 

xi>\i...>\l di>4...>di. (Al) 

A majorizes D if the following N conditions hold [2^ : 



Em 
fe=l 



Xt> 



^-^k=\ 



At 



di for m== l,...,iV- 1, (A2) 

Y..J^- (A3) 

In words: for each ttt, = 1, . . . , iV, the sum of the m largest 
elements of A is at least as large as the sum of the fc 
largest elements of D, with equality for m = N . 

Birkhoff showed that A majorizes D if and only if there 
is a double-stochastic matrix S such that 



D = SK 



N 



(A4) 



1 '^ij ~ 



where double-stochastic means that Sij > 0, ^ 

Let A be the sequence p3|) of initial eigenvalues. De- 
note by {\i)}1^Q~^ the set of eigenvectors for the initial 
Hamiltonian ^. Let D = {(«|ps+R(T)|i)}-i^"^ be the 
probability vector of energy level occupations in the final 
state ([5]). Eq. ^ implies (IA4|) with a double-stochastic 
matrix Sij = |(z|C/|j)|'', where the unitary operator U is 
defined by ^^. 

The Birkhof theorem implies that the sum of the 
largest k elements of D is not larger than that of A. For 



the purpose of cooling we want to make the sum of M 
largest elements of D as big as possible, and thus it has 
to be equal to the sum of M largest elements of A. This 
is realized if S permutes the M largest elements of A. 



Appendix B 

The example (|42|) was developed for M = 2n = 4. It 
can be generalized to arbitrary n: 

Po = TTo + 7ri + . .. + 7r„_i. A; = 0,. ..,n- 1, (Bl) 

T^2k = Po(7r/c + TTn+k), T^2k+1 = Pi (T^fc + 7r„+fe), (B2) 

W ^e{pi-pi) + ^^^^ ^ik{^Tk-^Tk)■ (B3) 



We assume 



Py^k EE /3/i > 1 



for 



k> 



Hence 



T^k>r 



0, 



Po -^ 1, 



(B4) 



(B5) 



EIL — J. r ^ ^IL— 1. 

, A^fc(7rfc - TTfc) H- >, , MfcTTfc (B6) 
k—l ^ — ^k—n 

En— 1 ^ — ^2n — 1 

, A*fe(^fc -""/c) + Ain > , T^k- (B7) 
k—l ^ — ^k—n 

Consider in detail the case n = 4, i.e., 8-level reservoir: 
Pi 



PW = piln 



■/3^i(pi7ro -TTi) 



1-Pi 



- [1 



- e-'^^i -t- e-'3''2 -I- e-''^^] 



[/3/ii(p 



1 -e 



-/3p 







-t- f3fi4{e-^^' + e-^^^)], (B8) 

where we used (jl21 I14p . We now minimize this expression 
over ^i < pL2 < /i3 assuming that /3/X4 ^ 1 is fixed [see 
(IB4[) ]. First we introduce new variables y2 and 2/3, 



/3^2 = hi 



PjH 

Z/2 



, ^/ia = In 



^M4 



2/3 



(B9) 



substitue them into (IB8I) and obtain after minimization: 



y2 = Poe 



'/3pi 



y3 = Pie 



-/Sa'i 



(BIO) 



Note that 2/2 > Vs seen from (jBlOp due to po > pi is 
consistent with /i2 < Ms- Next, we put (jBlOp back into 



PW = piln 



Pi 



+ 



I -Pi 

Puipi + e^^^'i ln[/3M4e^+'*[P«l 
1 + e-'^^i 



(Bll) 



where h[p] = —plnp— (1 — p)ln(l —p). Now the RHS of 
(IB11[) is to be minimized over ^i- This is done similarly 
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to dniHSll). Provided that lii[/3/X4ei+''[P°l] is sufficiently 
large, the maximization over /3/ii produced 



/3mi 



f3W = pi In 



Pi 



l-Pi 



In 



Pi In 



— ln[/3u4ei+''[P"l 
Pi 



Pi 



■ln[/3^4e^+''[P°l] 



(B12) 
,(B13) 



Hence the work needed for cooling scales double- 
logarithmically with maximal energy gap fj,^. The devi- 



ation of the ground-state probability (jBip from its max- 
imal value dSll) is controlled by dSSJ; see (|B9l [BTOI IBT2)) 
in this context. 



Eqs. (jB9[ IB12P show that the reservoir spectrum is 
self-similar: /3/i3 and /3/i2 depend logarithmically on /3/i4, 
while fii depend on /3/i4 doubly-logarithmically. 

Continuing this was one can show for a 2™"-'^ level 
reservoir the work can scale ln[. . . [ln[/3/i]] . . .] (logarithm 
is repeated m times) as a function of the gap. 



